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Abstract 

A very light scalar top (stop) superpartner is motivated by nat- 
uralness and electroweak baryogenesis. When the mass of the stop 
is less than the sum of the masses of the top quark and the lightest 
neutralino superpartner x?; as well as the of the masses of the light- 
est chargino and the bottom quark, the dominant decay channels of 
the stop will be three-body (3B: bW+Xi), four-body (4B: bffxi), 
or flavour violating (FV: ex?)- In this work, we investigate the direct 
and indirect constraints on a light stop, we compute the relative decay 
branching fractions to these channels, and we study the sensitivity of 
existing LHC searches to each of them. 



1 Introduction 



Supersymmetry can protect the scale of electroweak symmetry breaking 
from destabilising quantum corrections provided the masses of the Standard 
Model (SM) superpartners are not too large [T]. The most important states 
in this regard are those that couple the most strongly to the Higgs fields, 
namely the scalar top quarks (stops). To give the desired protection, the two 
stop mass eigenstates should be considerably lighter than a TeV [2]. 

A light stop can also play an important role in cosmology. The baryon 
asymmetry can be created within the minimal supersymmetric standard 
model (MSSM) by the process of electroweak baryogenesis (EWBG) if at 
least one of the stops is lighter than the top quark |3l HI |5]. A light stop 
can help to produce the observed dark matter relic density as well. When 
the lightest superpartner (LSP) is a mostly-Bino neutralino, the mechanism 
of thermal freeze-out tends to create too much dark matter. However, a 
light stop that is close in mass to the neutralino LSP can reduce its thermal 
abundance to the observed value by coannihilation [HI d El |9] . 

For these many reasons, an intense search for light stops is underway at 
the Large Hadron Collider (LHC), and a broad range of theoretical studies 
have been performed recently to find the most promising search channels pUf 
[m [HI [121 El Ha [m [13 HE] • The two decay channels that have been studied 
in the most detail are 

ti ^ tx?, h ^ bxt (2B) , (1) 

where Xi is the lightest neutralino and xt is the lightest chargino. When the 
stop is too light to decay in these channels, the dominant decay modes can 
involve three- or four-body final states such as [191 [201 [21] 

(3B,4B), (2) 

or have flavour violation [19l [20] , 

i^cxl (FV). (3) 

Of these, the two-body FV mode has been studied the most extensively [22l 
[23ll21l2Sl[2Sll2Zll2Hl[2nilSnil3Il[32llSSllMl[3S],andis typically assumed to 
be the exclusive decay channel when the three-body decay is kinematically 
forbidden. 
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In this work we study the LHC colhder signals of a very hght stop, with 
< rrit + itl^^q, in a comprehensive way. We place an emphasis on the 
scenario where only the 4B and FV decay channels are open. The branch- 
ing ratios of these modes depend on the underlying flavour structure of the 
theory, and we make the assumption that it is governed by Minimal Flavour 
Violation (MFV). In this case, we find that the 4B channel is frequently the 
dominant one despite its suppression by the large final-state phase space. To 
our knowledge, the LHC signals of a light stop decaying in this way have not 
been studied in detail. 

The outline of this paper is as follows. We begin in Sec. [2] by presenting 
the relevant stop interactions and the corresponding mass matrix, as well as 
our assumptions about MFV flavour mixing. Next, in Sec. [3] we investigate 
the indirect limits on very light stops from precision electroweak and flavour 
measurements. With these bounds in hand, we study the dominant decay 
channels of the light stop in Sec. |4j We then apply these results to make 
estimates of the sensitivity of existing LHC searches to a light stop. Finally, 
Sec. |6]is reserved for our conclusions. Some technical details can be found in 
the Appendix \K\ 



2 Stop Parameters 

To begin, we define our parameters and specify our assumptions about the 
flavour structure of the theory. In the absence of flavour mixing, the tree-level 
stop mass matrix is [1] 

M^-= ( "^Q3 + (i - lslr)c2^ml + ml mtX; 

with ttiq^ and m\jc the scalar soft squared masses, rrit mass of the top quark, 
tan/3 the ratio of Higgs vacuum expectation values (VEVs), and 

Xt = At- II* /i&Yi(3 (5) 

parametrizes the left-right (LR) stop mixing. A light stop mass eigenstate 
arises for small soft mass parameters or large LR mixing. 

Flavour mixing within the Yukawa couplings of the SM translates into 
some degree of flavour mixing in its supersymmetric extension. Additional 
sources of flavour mixing can arise from soft supersymmetry breaking pa- 
rameters, although these are strongly constrained by experiment [36l l37] . In 
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this work we will assume a specific flavour structure for the soft terms that is 
relatively consistent with existing flavour bounds, namely Minimal Flavour 
Violation (MFV) [38] ■ 

The structure of MFV is based on the observation that the SM has a 
global SU{3)^ flavour symmetry in the absence of Yukawa couplings. By 
promoting the Yukawas to spurious that transform under the global flavour 
group, the SM is made invariant under the symmetry. The flavour group is 
then broken only by the background values of the spurious. As applied to the- 
ories of new physics beyond the SM, MFV requires that this flavour structure 
continue to hold true, with any new sources of flavour mixing proportional 
to the Yukawa spurious. 

For the MSSM, the MFV hypothesis restricts the flavour structure of the 
soft parameters. The soft scalar mass-squareds become 6x6 matrices. To 
leading order in the Yukawa expansion they are given by [38| [39] 

ml = m\aiI + biYuY^ + b2YdYj + hYdYjY^Y^ + b,Y^Y^YdYj) (6) 



ml. = m\a2l + b,Y^Y^ + c^Y^Y^YiY^) (7) 

m% = Th\a^I + b(,YjYa) . (8) 

Similarly, the squark trilinear couplings are 

a, = A{a,l + brYaYj)Y^ (9) 

ad = Aia^l + bsYuY^)Yd . (10) 



The coefficients Oj, 6,, and q parametrize the expansion and are expected 
to be on the order of unity or less. Their values may be restricted within 
specific theories of flavour, but we treat them as free parameters at the 
superpartner mass scale |^ See Appendix |A] for a more complete description 
of our conventions for flavour and the soft masses. 

To compute flavour- violating stop decays, it is convenient to choose a 
speciflc basis for the Yukawa matrices where the up- Yukawas are diagonal, 

Yu = Xu Yd = VXd , (11) 

where V is the CKM matrix, and \u,d are the diagonal up- and down- Yukawa 
matrices. 

"'^This is consistent with but more general than the commonly- used assumption that 
the flavour mixing vanishes at some high-energy scale. 
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For our study we choose a set of fiducial soft parameters that allows us 
to isolate a light stop state. We take rh = 1500 GeV and ai = 02 = as = 1, 
which has the effect of decoupling the first- and second-generation squarks, 
and we choose similar soft masses for the sleptons. Much smaller values of 
(^q)33 ('^ciOas obtained by suitable choices of bi and 65, and these 
translate into one or more light squark states. Since these states only contain 
very small admixtures of the first- and second-generation squarks, we identify 
them with a pair of stops and a (mainly left-handed) sbottom. The rest of 
the bi and q parameters are scanned over the range [—1, 1]. We also set the 
gluino soft mass to M3 = 1000 GeV, but we consider a range of values for 
the electroweak gaugino soft masses Mi and M2 and the Higgsino mass fi. In 
implementing these values within the spectrum generators SuSpect 2.41 [ID] 
and SoftSUSY 3.3.5 [3T]|, we input these parameters as running values defined 
at Ms = y/mtrriQ. Finally, we consider the specific representative cases of 
tan^ = 10, 30. 

3 Limits on a Light Stop 

Stops have been searched for extensively at the Tevatron and the LHC A 
light stop can also influence a number of SM observables by its contribu- 
tions in loops. We review here the existing limits on stops from direct col- 
lider searches and we study the indirect limits on a light stop from precision 
electroweak measurements and flavour violation. We also comment on the 
constraints imposed by the rates of Higgs production and decay. 

3.1 Direct Collider Searches 

The ATLAS and CMS collaborations at the CERN LHC have succeeded in 
confirming the SM to such an extent that squarks of the first and second 
generations are constrained to be heavier that a TeV H3] . More recently, 
these collaborations have turned their attention to investigating the squarks 
of the third generation created both directly and indirectly. Recent studies 
rule out direct stop pair production for stop masses up to nearly 600 GeV 
under the assumption that the only decay is t — )■ txi with m^o — > and 
Xi a stable LSP |S1 US]. The limits weaken for larger values of the LSP 
mass, but remain in the range of many hundreds of GeV. Searches for direct 
sbottom pair production with b — )■ 6x1 yield similar bounds on the lightest 
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sbottom mass [IS]. Related searches for direct stop production with the 
decay i bxt together with xt ~^ W'^^*^Xi also give exclusions as large as 
mi > 450 GeV [HI HS]. Finally, light stops or sbottoms can be created in 
gluino decays. Assuming one or the other of the topologies g — )■ ttxi with 
i — )■ txi or ^ — 7- 66 with b — )■ bxi, the constraint on the gluino mass extends 
up to m-g = 1.3 TeV [SI SB]. 

Collider limits on stops and sbottoms can be much weaker when they 
decay in other ways, such as the 3B, 4B, and FV channels. Searches at 
LEP rule out a stop decaying in the FV channel up to about 100 GeV [19], 
and we expect similar limits to hold for the 3B and 4B decays. At the 
Tevatron, specific searches for a light stop in the FV mode give limits up to 
mi ^ 180 GeV, but these limits disappear when the stop becomes moderately 
degenerate with the LSP [501 [51]. Almost no limit is expected when the 3B or 
4B channels are dominant [23]. The LHC collaborations have not performed 
dedicated searches for a stop in the 3B, 4B, and FV modes. The goal of 
this work is deduce the sensitivity of existing LHC searches to these stop 
channels. 



3.2 Indirect Bounds from Electroweak and Flavour 

The most important effect of a light stop on precision electroweak observables 
is to modify the prediction for the mass of the W boson, corresponding to 
a shift in Ap (or T in the STU parametrization of oblique parameters [52]). 
This effect has been studied recently in Refs. [531 Ell ES] in light of the 
improved determination of Mw [56] together with evidence for a Higgs boson 
of mass mh = 125 GeV [SZl EH]- These results imply Ap = (4.2 ± 2.7) x 

10-^ [S1E5]. 

The change in Ap due to light stops is given by [59] 



Ap ^ [-^?c'^o(m? ,m|) + clFoimlml^) + s?Fo(m|,m?j] , (12) 

where 9i is the stop mixing angle and Fo{x,y) is a loop function, 

Fo{x,y)=x + y-^^\n(-] . (13) 



x-y \y 

We apply this formula in the present case where we have squark flavour 
mixing dictated by the MFV paradigm. Since the inter-generational mixing 
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is numerically small, we are able to identify specific squark mass eigenstates 
with ti, t2, and Bl, and apply this formula directly. Additional corrections 
due to flavour mixing are expected to be highly subleading. We also neglect 
the contributions from neutralinos and charginos, which are also expected to 
be less important [HI EHl EI]- 

A light stop can also contribute to flavour mixing, and to the branch- 
ing fraction BR(i? — Xg'y) in particular. This effect is present even if 
the scalar soft parameters are all diagonal in the super-CKM basis (dis- 
cussed in Appendix |A]) due to loops involving a top and a charged Higgs or a 
stop and a chargino. These new contributions vanish in the supersymmetric 
limit, but can modify the branching fraction when supersymmetry breaking 
is present [621 EHl El]. In our analysis, we use the HFAG experimental re- 
sult BR{B X,7) = (3.55 ± 0.24 ± 0.09) x 10"^ [65] together with the SM 
prediction of BR{B X,-f) = (3.15 ± 0.23) x 10"^ from Ref. [M]. We also 
consider limits from BR(5, n+i^-) [67| and Ao{B K*-f) [33 |65]. 

We investigate the contributions of a light stop to BR.{B — )■ ^^7) and 
other flavour observables by using SoftSUSY 3.3.5 [41j to compute the super- 
partner mass spectrum and mixings and passing the output to SuperlSO 3.3 [68] 
to evaluate the modifications to flavour observables. Input and output are 
achieved in the SLHA2 format [69], and a detailed description of the interface 
is given in Appendix |A| 

With these tools, we have performed a series of scans over stop param- 
eters by varying the A, bi, and q MFV parameters, while holding fixed 
= 1 and = (1500 GeV)^ to maintain heavy first- and second-generation 
sfermions near 1500 GeV. We have also set M3 = 1000 GeV, Mi = 200 GeV, 
M2 = 350 GeV, and fi = Ma = 500 GeV. The resulting allowed regions, 
consistent with both precision electroweak and flavour constraints at the 
2(j-level|^ are shown in Fig. fl] for different values of the stop mass m^-^, 
the stop mixing angle | cos 9f:\ FT and the stop mass-matrix mixing parameter 
Xt = Xt33 . In both panels, we have fixed tan /3 = 10 and we have restricted 
100 GeV < m^^ < 300 GeV and m^^ > 100 GeV. We also show the ef- 
fects of imposing additional restrictions on the mass of the lightest sbottom: 
mi_^ = 300 GeV, 600 GeV, which are motivated by the strong limits placed 
on sbottoms by the recent analysis of Ref. [46J . 

Our primary conclusion from Fig.[T]is that electroweak and flavour bounds 

^The ranges are Ap e [-1.2,9.4] x 10"'' and BR{B X,-f) g [2.77,4.33] x 10"^. 
■^Our convention is \cf\ — 1 (0) when the hghtest stop is purely left (right) -handed. 
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Figure 1: Allowed light stop mass and mixing angles for tan/3 = 10. 



alone permit a broad range of mixing angles for a light stop, but that de- 
manding a heavier sbottom forces the light stop to be mainly right-handed]^ 
Despite the broad range of allowed mixing angles, we also find a severe re- 
striction on the stop mixing parameter Xt for smaller values of '^'q^^, which 
comes partly from Ap but mainly from BR(i? — )■ Xsj). We have also ex- 
amined the allowed ranges for tan/3 = 30, and we obtain similar results, 
albeit with a stronger restriction on Xt. The limitations on Xt derive from 
an interplay between the chargino-stop and charged Higgs-top contributions 
to BR(i? — 7- Xg'-f). These must cancel to some extent when there is a light 
stop in the spectrum, with the specific condition depending mainly on /i. Ma, 
tan/?, and Xt. Even so, at the expense of some amount of tuning between 
these parameters, a broad range of stop properties is possible. 

Let us also mention that among the bounds from flavour mixing, we 
find that BR.{B — )■ Xg'y) gives the most significant limit for the param- 
eters we studied, and is more restrictive than BR(Ss — )■ |SZ] and 
Ao(-B — )■ K*'y) [65]. We also find that the contributions to flavour observables 
are typically dominated by the CKM mixing already present in the super- 
CKM basis. More precisely, the deviations from flavour-diagonal squark mass 
matrices in the super-CKM basis allowed by MFV do not appreciably change 
the effects of the superpartners on the flavour observables we have studied. 
Indeed, in the case of BR(i? — )• Xs'y), we find substantial agreement with 
SuSpect 2.41 |10] which does not include any additional (non-super-CKM) 
flavour mixing. 

A light sbottom might still be consistent with LHC searches if it is highly degenerate 
with both the light stop and the neutralino LSP 05]. 
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3.3 Higgs Production and Decay Rates 

Light stops are also closely linked to the mass of the SM-like Higgs boson 
in the MSSM. For many of the sets of stop masses we consider in this work, 
the mass of lies well below the mass of the tentative Higgs resonance near 
125 GeV [701 EI]- Within the MSSM, this discrepancy can be resolved by 
pushing the left-handed stop soft parameter ^^^33 to values well above the 
TeV scale [31 172] . This problem can also be addressed by expanding the field 
content to include an additional singlet, as in the NMSSM. In both cases, 
these fixes can be achieved without necessarily altering the collider signatures 
of the light stop. 

A light stop can also modify the production and decay rates of the /i° 
Higgs [73]. When the light state is mainly unmixed, the effect is to enhance 
the production rate in gluon fusion and to suppress the branching fraction 
to photon pairs [711 ESI ESJ. Conversely, when the the light stop is highly 
mixed it has the opposite effect, pushing down the rate for gluon fusion and 
increasing the branching ratio to diphotons [77]. The Higgs rates can also 
be modified by the presence of other MSSM (or beyond) superpartners that 
run in the loops [721 EH] or that mix with the Higgs [791 EO]- Again, these 
extensions of the MSSM need not modify the collider signals of the light stop. 

Light stops can also contribute in a more direct way to Higgs searches 
by the formation of stoponium, a stop-antistop bound state [HD EH ESI El] ■ 
Stoponium decays can mimic those of a Higgs boson. However, for m^^ > 
m^o, the dominant stoponium decay mode is frequently to a pair of Higgs 
bosons (particularly if the stop state is not highly mixed). Current LHC 
searches do not appear to be sensitive to this decay channel [83]. 

For these various reasons, we do not constrain the spectrum of light stops 
based on searches for and the tentative discovery of the Higgs boson. 
Instead, we focus on the direct LHC limits on a hght stop state. 

4 Light Stop Decay Modes 

We consider the decays of a light stop to the lightest neutralino or chargino 
in the spectrum, taking the rest of the superpartners to be heavy. We also 
focus on the case where t txi is kinematically forbidden. The stops then 
have four possible decay modes: 

1. two-body (2B): bxt ■ 
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2. three-body (3B): bW+Xv 



3. four-body (4B): bf/'xl- 

4. flavour- violating (FV): cXi- 

Our analysis shows that each of these modes tends to dominate in a different 
kinematic region, depending on the stop-neutralino and stop-chargino mass 
differences. 



4.1 Flavour- Violating Stop Decays in MFV 

The most model-dependent of these decays modes is the FV channel. We 
apply the results of Ref. [85] to calculate the FV decay width within the 
framework of Minimal Flavour Violation (MFV) [38] • The general form of 
the stop-charm-neutralino coupling is [85] 

-^cixS = (^(VlPr + yRPL)x\i + h.c. , (14) 
in terms of which the FV decay width is given by 

where Y = v^|i/lP + Ii/rP- 

The couplings and can be computed reliably in the mass-insertion 
approximation. The relevant flavour mixings in this case arise from the 
following off-diagonal elements of the soft masses [85] : 

(m|j23 = m%{YdYj)23=rh%2XlVcbV,l 

irhlj23 = m^ci{Y^YdYjY^)23 = m^ciXtXcXlV,tV;, (16) 

(a„)23 = Abr{YdYjY^)2s = AbrXtXlV^tV;, 

(a„)32 = Ab,iY,YjY^),2 = AbrXAlVcbV;,. 

In writing these expressions we have neglected terms proportional to the 
charm Yukawa. These contributions are negligible unless the LR stop mixing 
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is tuned to nearly zero. With these off-diagonal soft terms, the stop-charm- 
neutralino couplings are 



VL = ^i^^'^i {rh%c, + AvJjX,Si){-/-Nn - -^N,,), (17) 

VR = 0(Ac) (18) 

where Cf = cos 6j: and sj; = sin 6^ correspond to the LR stop mixing, Vu is the 
up-type Higgs expectation value, and is the neutralino mixing matrix. 
Note that the relevant contributions to the FV decay width correspond to 
the stop mixing into the left-handed charm squark c^. All the contributions 
involving the right-handed Cr squark are proportional to the charm Yukawa, 
which are subleading and are neglected here. 

We compute the FV decay widths by calculating the mass spectrum (with- 
out flavour mixing) with SuSpect 2.41 ^40J. Comparing with SoftSusy 3.3.5, 
corrections to the masses and the LR mixing from off-diagonal soft terms are 
small provided the flavour-diagonal terms match up. Using these masses, we 
then apply Eqs. ( [I5|l7[18[ ) to obtain the decay width. 



4.2 Other Decay Modes 

For the 2B, 3B, and 4B modes, we use SuSpect 2.41 |10] to calculate the 
mass spectrum and we interface it with SUSY-HIT [86] to find the decay 
widths. While the 2B and 3B channels are treated in full, the SUSY-HIT 
implementation of the 4B mode is an approximation to the full result [21] 
where the outgoing //' fermion decay products are treated as massless. We 
expect this to be a good approximation except for very close to the kinematic 
boundary where we apply a phase space factor to account for the fermion 
masses. 



4.3 Phase Diagrams 

With the decay widths in hand, we can chart out the dominant decay channels 
in different regions of the underlying parameter space. We find that the most 
likely decay channels are determined largely by phase-space considerations, 
which depend in turn on the stop, chargino, and neutralino masses. When 
the 2B and 3B channels are disallowed, we also find that the 4B mode can 
overcome the FV mode within the MFV framework. 
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Figure 2: Dominant stop decay modes of a stop with mass = 225 GeV for 
tan/3 = 10, with cos^^ = 0.1 (left) and cos^j = 0.7 (right). The four regions 
correspond to dominant decays by the 2B (gray), 3B (red), 4B (yellow) and 
FV (green) channels. We also show the kinematic boundaries between these 
decay channels by solid black lines. 



In Fig. |2] we show the dominant decay channel of the lightest stop in the 
Mi—fi plane for a stop mass of = 225 GeV, tan /3 = 10, and mixing angles 
of cosOf = 0.1 (left) and cos6^ = 0.7 (right). The four regions in the plots 
correspond to dominant decays by the 2B (gray), 3B (red), 4B (yellow) and 
FV (green) channels. The phase space boundaries between the various decay 
channels are shown with solid black lines. In these plots we have also fixed 
Ma = 350 GeV, Ma = 500 GeV and 62 = &7 = 0.5. All other input terms 
are as described at the end of Sec. M 

From these plots we see that the 2B channel dominates for smaller values 
of /i, when mi > ruh + and the decay is open. This is to be expected, 
since the channel is neither suppressed by phase space nor flavour mixing. 
As this mode turns off, the 3B channel comes to dominate at smaller values 
of Ml where it is allowed, up to > m?, + mw + ^x?- -^^^ larger Mi, the 
3B channel shuts off as well and the AB and FV channels take over. The 
boundaries between these regions are also nearly identical to the phase space 
boundaries. Note that the AB channel is cut off when (m^ — m^o) < rrih, in 
which case the FV mode is all that is left. 

The interplay between the 4B and FV decay modes is more complicated. 
While the 4B mode has significant suppression from the four-body final-state 
phase space, the FV channel is reduced by the level of flavour mixing. Very 
approximately, the relative suppression factors are 1/4! (47r)^ ~ 2 x 10~^ for 
the 4B decay and (mfetan/3/t;)^|K6p ~ 5 x 10-6(tan/3/10)l However, FV 
decays can receive further suppression if the light stop is mostly right-handed, 
or they can be enhanced for larger values of tan (3. 
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Figure 3: Dominant stop decay modes of a stop with mass = 225 GeV for 
tan/3 = 30, with cosO^ = 0.1 (left) and cos9^ = 0.7 (right). The four regions 
correspond to dominant decays by the 2B (gray), 3B (red), 4B (yellow) and 
FV (green) channels. We also show the kinematic boundaries between these 
decay channels by solid black lines. 



Comparing the left and right panels of Fig. [2| the enhancement of the FV 
channel relative to 4B with increasing component is evident. When the 
lightest state is mostly right-handed, as in the left panel of the figure where 
cos6^ = 0.1, the 4B decay dominates over the FV nearly up to the kinematic 
limit. For this reason, together with the mass hierarchy between the stop and 
the scharm, we obtain a larger branching fraction in the 4B mode relative to 
previous analyses [201 IHI| ■ 

When the stops are maximally mixed, cos^^ = 0.7, the FV channel takes 
over well beyond the kinematic boundary. This trend continues as the left- 
handed content of the light stop is increased further. However, this also gives 
rise to a light sbottom squark and a new decay channel i — )■ bW~^^*^ that can 
dominate in what would otherwise be the 4B and 3B regions. Given the very 
strong direct-search constraints on a light sbottom as well as the indirect 
bounds from precision electroweak and flavour mixing, we do not consider 
this possibility here. 

The FV mode can also be enhanced relative to the 4B mode for larger 
values of tan /3. In Fig. [3] we show the dominant stop decay channels for 
nearly the same parameters as in Fig. [2| = 225 GeV, M2 = 350 GeV, 
b2 = h'j = 0.5, but now with tan/3 = 30 (instead of tan/3 = 10). In the left 
panel we have cosOi = 0.1, and cos% = 0.7 in the right. As suggested by 
our simple estimates, the FV channel dominates over a significantly larger 
region than before. We conclude that, in MFV, either the FV or 4B decay 
modes can dominate in this region of the superpartner spectrum, and it is 
important to search for both. 

The absolute decay widths obtained in the FV-4B region receive a sig- 
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nificant amount of suppression. Even so, we find that they are generally 
prompt on collider time scales unless the stop-LSP mass difference becomes 
very small, in agreement with Ref. [85]. For our collider studies, we will as- 
sume that all the decays are prompt. However, the emergence of a displaced 
stop decay vertex could allow for the direct measurement of the underlying 
flavour structure and would be interesting to study in its own right. 
On the other hand, a stop that is long-lived on the time scale of the LHC 
detectors is firmly ruled out to masses well over m^^ > 500 GeV [88] . 

5 LHC Search Sensitivity 

In this section we investigate the sensitivity of existing LHC searches to 
a very light stop decaying via the FV, 4B, or 3B modes. We make use of 
public ATLAS and CMS analyses to estimate exclusion limits in the m^—m^o 
plane for stop masses in the range 100-350 GeV. Although these experimental 
studies were not optimized for light stops, we find that they imply quite 
stringent bounds on such a state. 

Throughout our investigation, we treat the light stop within a simpli- 
fied model where it decays exclusively (and promptly) via the FV, 4B, or 
3B modes to a neutralino LSP and SM particles within the corresponding 
kinematically relevant region of the — m^o plane. As seen in the previous 
section, varying the LR stop mixing, tan/3, or the amount of flavour mix- 
ing can modify the relative branching ratios of the FV and 4B channels]^ 
Exclusion limits for intermediate cases can be obtained by interpolating be- 
tween the separate FV and 4B cases studied here. We do not consider light 
sbottoms, or production via gluinos or other superpartners. 

5.1 Event Generation and Simulation 

We generate pp — > tt* events at ^/s = 7 TeV using MadGraph5 vl.4.3 [5U] . 
The parton-level output is passed to Pythyia 6.4 |90j for decays, showering, 
hadronization, and matching of parton showers to matrix elements. For 

^ While these decay modes are usually prompt, they are slow enough that the stop will 
hadronize before it decays. We do not attempt to account for this effect, but we expect 
that it will only be significant close to the kinematic cutoffs based on the typical momenta 
of the decay products. 



14 



matched samples, the MLM scheme is used with shower-A;_L and Qmatch = 
50 GeV. 

Simulation of the LHC detectors is performed with Delphes [91], using 
the anti-fcy jet reconstruction algorithm. ATLAS and CMS detector/trigger 
cards are applied to obtain a more realistic comparison of our analyses to 
the respective LHC searches. We have cross-checked our results with PGS4 
(using cone jet reconstruction with A R=0.7) |92]. The differences in the 
distributions of hardest jet pt and missing energy {^t) are within generally 
within 5-10% for the two detector simulators. In terms of acceptances to 
various searches, Delphes tends to be somewhat lower than PCS. 

In generating stop samples for our analysis, we are faced with a practical 
challenge. Many of the existing LHC SUSY searches are geared towards 
finding relatively heavy coloured states that decay to a much lighter LSP. 
These searches frequently require multiple hard jets and apply strong cuts 
on or riieff = T,jPT + ^t- such, they are typically only sensitive to the 
FV or 4B stop decay modes in certain specific regions of the stop-production 
phase space. 

Since we are investigating relatively light stops, the typical energy scale 
in the events tends to be somewhat small. Furthermore, the momenta of the 
stop decay products are distributed between multiple final states in the 4B 
channel, and are also limited in the FV channel when the stop-LSP mass 
difference is small. Thus, the light stop events usually produce soft decay 
products and small amounts of ^j.^ particularly if the LSPs are emitted back- 
to-back. For these light stop events to pass many of the LHC SUSY analysis 
cuts, they must be produced in conjunction with one or more hard jets. In 
addition to increasing the total energy in the event, the radiated jet gives a 
net boost to the tt* system which can lead to an observationally significant 
amount of missing energy [93l IM] . 

The challenge, therefore, is to simulate efficiently this small corner of 
the stop production phase space where one or more hard jets are emitted 
without being swamped by the much larger portion where the additional 
jets are relatively soft. For this, we investigated the following generator-level 
samples: 

tt* , tt* + Ij (unmatched), tt* + Ij (matched), 

with matching scale Qmatch = 50 GeV for Ij matched, and requiring pt > 
50, 100, 150, 200, 250, 300 GeV for the hardest jet at the parton level in the 
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unmatched samples. Of these, the matched sample is expected to best repro- 
duce the full pt distribution of the additional jet [95]. However, each given 
sample is strongly dominated by events where the extra jet lies near the 
matching scale, leading to poor statistics in the higher-p-p region of greatest 
interest. 

To improve the statistics of the samples we generate, we compared the 
matched sample to a set of unmatched tt* + Ij samples constrained to have 
a jet Pt greater than 50 — 300 GeV. For our comparison, we studied the 
relative signals predicted by these samples for a range of fiducial kinematic 
cuts typical of the ATLAS and CMS analyses considered in this paper: 

> 130, 200, 300 GeV m^// > 500, 700, 1000 GeV 

Ht > 350, 500, 600 GeV prji > 80, 130 GeV . 

We find that the tt* + Ij unmatched sample with prji > 150 GeV provides 
a good balance between obtaining useful event statistics {i.e. a reasonable 
fraction of the generated events pass the kinematic cuts) and not artificially 
eliminating the signals. For example, the unmatched sample with Prji > 
100 GeV leads to the same effective signals to within 10 — 20%. 

We have also compared the pt distribution of hardest-jet in the Ij un- 
matched sample with prji > 150 GeV to matched samples with one or two 
additional hard jets. Relative to the matched Ij sample with Qcut = 50 GeV, 
we find excellent agreement for Prji > 150 GeV (as is expected). Compar- 
ing to the matched 2j sample (matched with Ht > 100 GeV), we also get 
excellent agreement for the Ptji distributions. While the 2j matched sam- 
ple tends to produce a harder second jet, we find that including it does not 
appreciably increase the effective signals that pass our sample kinematic cuts. 

The stop production sample we generate arises at NLO and represents 
only a small subset of the total stop production phase space. As a result, 
we do not attempt to apply a K-factor to the cross section obtained by 
comparing the inclusive LO and NLO cross sections (which is about 1.4 in 
this case based on Prospino2.1 [96\). Instead, we estimate the effect of the 
theory uncertainty in the cross section on the LHC exclusions by simply 
varying the MadGraph production cross sections by ±50%. 

5.2 LHC Searches and Exclusions 

The LHC analyses that we apply to our light stop samples are summarized 
in Tables fTl and [21 These searches constrain different regions of the — m^o 
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ATLAS Analysis iC/fh-^) 


Final states 


Decay channel 


Mnnnipt {\ DDI PTtI 


< l-S? 0/ 


FV 4R 


Tpfoiir^ C\A\ no] 


> 2-47 0/ 


FV 4B 

J. V ^ ^ 1 J 


Monojet (4.7) [98] 


<2j,0/ 


FV, 4B 


j + b + :^T (0-83) [99] 


> 3j, > 16, 0/ 




j + l + ^T (1-04) [lOOj 


>3-4j, 11 




j + b + l + ^T (1-03) [lUU 


>4j, > 16,1/ 




j + b + l + ^T (2.05) [102J 


> 3-4j > 16, 0-1/ 




^ // (4.7) [ing 


> Ij, 1^ 




J + / + 6 + (4.7) [M] 


>4j,16,l/ 




j + 6 + ^'r(4.7) 1105J 


> 4j, > 26, 11 





Table 1: Summary of ATLAS searches at ^/s = 7TeV that constrain the FV, 
4B and 3B decay modes. 

plane, and taken together, provide significant bounds on a light stop. All the 
searches considered have ^/s = 7 TeV. We defer an analysis of the recent 
SUSY search results obtained at y/s = 8 TeV, which appeared as this paper 
was nearing completion, to a future study. 

We consider stop decays in the FV, 4B, and 3B modes, and the limits 
from each LHC analysis are applied separately. In every case, we impose the 
selection requirements described in the corresponding report to our simulated 
signal and we compare to the listed backgrounds (accounting for the stated 
background uncertainties) to obtain limits on the stop and neutralino masses. 
The corresponding 95%c.l. exclusions are summarized in Fig. |4] (FV, 4B) 
and Fig. [s] (3B). In these figures, we show the result for the central value 
of the MadGraph-generated cross section with solid lines, and the exclusions 
when this cross section is varied by ±50 % with dashed and dotted lines. 

Note that many of the searches listed in Tables [T] and [2] do not yield 
a constraint. This fact is shown in the rightmost column of the Tables, 
where we list the decay channels for which the corresponding searches are 
relevant; when no channel is listed, no interesting bound was obtained. In 
the discussion to follow, we will focus on the LHC analyses that do lead to 
limits. 
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CMS Analysis (C/fb-^) 


Final states 


Decay channel 


Mnnnipf (11) ll Dfi 


< 9i 0/ 


FV 4R 

X V 5 TiXJ 


7 + Err, (1 1) 11071 


> 3? 0/ 




fyrr fl 141 11081 


> 2i Q-lKii) 0-1^ 






> 3-47 > Iftl/ > 1/) 


FV 4B 

X V ^ ^ 1 J 


Razor (4 4) IllOl 


> 2? 0-2/ 


FV 4B 3B 


J+/ + ^'rfl.l) [111] 


> 3-4j, 1/ 




(1.1) 11121 


> 3-4j, > 0-16, 0/ 


FV, 4B, 3B 


Monojet (5.0) [113] 


<2j,0/ 


FV, 4B 


J + ^^(4.98) 113] 


>3j,0/ 


FV, 4B 


ar (4.98) pjL4j 


> 2j, 0-36, 0/, 




j + 6 + (4.98) [115J 


> 3j, 1-36, 0/ 


FV, 4B, 3B 


Razor (4.7) [116J 


> 2j, > 16, 0-2/ 


FV, 4B, 3B 


Mt2 (4.73) [117J 


> 3-4j, > 0-16, 0/ 


FV, 4B, 3B 



Table 2: Summary of CMS searches at ^/s = 7TeV that constrain the FV, 
4B and 3B decay modes. 

5.2.1 FV mode 

The FV mode has t — ?■ cXi, which suggests that it will be probed most 
effectively by searches for jets and missing energy. Thus, we apply the ATLAS 
and CMS monojet, jets-f-^-^. Razor, and Mt2 analyses to it. However, due to 
a sizable mistag rate for c jets to be identified as 6 jets (about 10 %), we find 
that the FV mode is also constrained by several 6-jet analyses. In contrast, 
we find that searches with leptonic final states are generally not sensitive to 
this decay channel. The regions of the stop-neutralino mass plane excluded 
by the analyses we consider for the FV decay mode are shown in the left set 
of panels in Fig. |4j 

Monojet searches require at least one hard jet and missing energy, and 
they are good probes of the small — m^o region where the outgoing charm 
quark is soft and difficult to detect [22]. In this case, the hard jet arises 
from QCD radiation. ATLAS and CMS monojet searches apply various sets 
of cuts on the transverse momentum of the hardest jet ptj^ and the missing 
energy ^j.. They also veto on events with hard leptons, as well as events with 
too many additional hard jets [971 EHl I106[ 1113] . We find that the harder cuts 
on ptji and used by these analyses, which range up to 500 GeV, tend to 
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be the most effective at probing the stop signaL Among the set of events that 
pass the selection requirements, the stops must be created with a very hard 
extra jet. This can give the stop pair a large pr, which leads to significant 
missing energy but can also make the charm jets visible. From this point 
of view, the less-stringent veto on additional jet activity in Ref. [Ml 1113] 
relative to Ref. [HZ] is beneficial]^ 

The presence of one or more hard QCD jets and the resulting boost of 
the stops can lead to multiple hard jets in FV stop events. Thus, searches 
for multiple jets and missing energy can also be sensitive to light stops. The 
recent ATLAS jets+^j" analysis of Ref. |12] divides their search into 2 — 4j 
channels, and imposes a lepton veto and cuts on and rrieff = YlijPT, + 
^j.. The rrieff and requirements of this analysis are geared towards 
heavier superpartners and are quite severe, with me// typically greater than 
1000 GeV, limiting their sensitivity. In fact, the more recent ATLAS analysis 
of Ref. |104j with 4.7 fh~^ of data imposes even stronger requirements on m-e// 
and is less sensitive than Ref. [12], despite having much more data. The 
kinematic requirements for the corresponding CMS searches for jets-t-^'^ are 
less severe, and give slightly stronger exclusions. The CMS search for jets and 
missing energy of Ref. [33] requires at least 3 jets, Ht = YlijVTj > 500 GeV, 
missing energy, and no hard leptons. This search places strong limits on 
cross-sections for new physics, excluding a region of the — m^o plane that 
is complementary to the monojet analyses. 

The FV mode is also constrained by searches involving b jets due to the 
significant (~ 10%) mistag rate of c-jets. The minimal requirements of the 
CMS jets-t-b-f-^-j. search |109] . at least three jets with px > 50 GeV and a 
b tag, can be satisfied in this casej^ This analysis also has weaker cuts on 
Ht and missing energy (limited via the variable) relative to the jets-f-^^^ 
searches, which leads to a significant exclusion of the FV mode. Let us also 
point out that the observation of a light squark in multiple light- and heavy- 
quark flavour channels could provide information about the flavour structure 
of the stop decay. Even though the c quark is misidentified in this case, the 
decay could still be differentiated from other flavour-violating channels such 
as t ^ uxi, (which is negligible in MFV), since the mistag rate of w-jets as 
6-jets is about 1%. 

^Thc bounds from the monojet searches are also Hkely to be vahd when the stop decay 
is sUghtly displaced. 

^ A small fraction (about 5%) of the stop events passing the cuts have genuine b jets 
produced in the parton shower. 
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Light stops are also constrained by the CMS Razor and Mt2 searches. 
The Razor variables provide a way to use event-shape information to reduce 
backgrounds while not imposing overly strict kinematic requirements |110j, 
I116j . The stop signal populates mainly the Mr < 1000 GeV region, and the 
analysis of Ref. |110] restricts a significant part of the stop-LSP mass plane. 
The CMS M'j'2 analyses of Refs. \11'2\ I117j have two separate signal regions, 
with and without a 6-tag requirement. Applying the analysis of Ref. [112] . 
the HighMT2 search channel (> 3j, 06) excludes a greater region than the 
LowMt2 channel (> 4j, > 16) due to the 6-tag requirement of the latter. 
Similar exclusions are found for the more recent search of Ref. |117] . While 
this analysis has considerably more data than Ref. |112j . it also imposes more 
stringent kinematic requirements leading to about the same exclusion. 

Our results for the FV mode appear to be compatible with previous anal- 
yses of this channel in the monojet and jets+^j. channels |2H |25l |26l |27] . 
although Ref. [33] obtains a somewhat stronger limit from the jets-f^^" chan- 
nel than we do. The limits obtained in our analysis from jets+b+^^y and 
MT2/Razor analyses provide new and slightly stronger exclusions. 
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Figure 4: ATLAS and CMS exclusions on FV mode (left) and 4B 
mode (right). The two straight, diagonal lines correspond to = m^o (up- 
per, FV), rrif = m^o + rrih (upper, 4B), and = m^o + nib + m\Y (lower, 
both). The upper lines are the kinematic limits for these modes, while below 
the lower line we expect the 3B channel to take over. 
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5.2.2 4B mode 



Each 4B stop decay produces the final state t — )■ bf f'xi, where the //' 
fermions come from the decay of an off-shell W^. Events with a pair of stops 
can therefore produce zero, one, or two leptons. While the potential presence 
of leptons might seem promising, we find that they tend to be somewhat 
soft and difficult to identify when produced by this decay modej^ Hadronic 
analyses turn out to be more effective, although their reach is slightly weaker 
than for the FV decay channel. 

Among the primarily hadronic search channels, we find that the ATLAS 
and CMS monojet and iets+^j^ analyses are moderately effective. The mono- 
jet searches are sensitive mainly to the case where the hard jet is produced 
by QCD radiation, and as a result the 4B region excluded by the monojet 
searches is similar to the FV case. The ATLAS jets+^j^ search of Ref. [12] 
requiring two hard jets is effective for probing the small (m^- — m^o) region. 
However, the higher number of final states in the 4B decay results in softer 
jets relative to FV decays and a smaller exclusion region. We find a similar 
result for the CMS jets+^rp search of Ref. [43j . 

The CMS Razor analyses |110[ 1116] constrain the hadronic 4B mode only 
in the hadronic (HAD) channel. Although some of the Razor analyses also 
search for one or two e or in the final state, the leptonic 4B channel is not 
picked up by the search. The light stop decays are constrained by the low 
transverse mass search region, Mr < 1000 GeV, and as a result, the excluded 
region is similar to the FV mode. The more recent search of Ref. |116] has 
an additional 6-tag requirement and more data, and excludes a larger region. 
Similarly, the CMS Mt2 analyses of Refs. |112l 1117] that demand a 6-tag 
are also effective. However, the net exclusions of these analyses are slightly 
weaker than for the FV mode, due to the higher multiplicity and softer 
momenta of the 4B stop decay products. Searches for 6-jets and missing 
energy at CMS (with a lepton veto) are also effective at probing the 4B 
decay piMITT^ . 

5.2.3 3B mode 

The 3B mode is sensitive to analyses with &-jet final states, in particular 
the CMS 6-jet, Razor, Mt2 analyses. The newer CMS 6-jet search |115] 
constrains the 3B mode, with an exclusion region comparable to the FV and 

^Thc requirement of multiple hard jets in these analyses is also problematic. 
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m^(GeV) m^(GeV) 




Figure 5: CMS exclusions on 3B mode. The upper solid diagonal line shows 
the kinematic limit for this decay, = m^o + mb + mw, while the lower solid 
diagonal line has m^- = m^o + nit where the 2B decay channel is expected to 
take over. 



4B modes. In particular, the signal region with one loose &-jet, IBL, is most 
sensitive to the decay mode. 

The CMS Razor analyses |110[ I116j are sensitive to the 3B mode in the 
hadronic search channel (HAD), with stronger limits from the more recent 
search. The decays are constrained by the transverse mass search region, 
Mr < 1000 GeV. The CMS Mt2 analyses [HSl [HI] provide the most strin- 
gent constraints, up to ~ 300 GeV. In particular, the 6-tag signal region 
is most sensitive to the 3B decay mode. 
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6 Conclusion 



Very light stops can be consistent with current experimental bounds when 
they decay in ways that are challenging to detect at the LHC. In this study 
we have attempted to investigate the properties of light stops in a system- 
atic way, and to estimate the bounds placed on them by indirect precision 
measurements and existing searches at the LHC. 

When a stop is very light, it can decay according to i ^ cXi (FV) and 
t — )■ bff'xi (4B). The relative branching fractions depend on the gauge 
content of the light stop, the value of tan /3, and the amount of flavour mixing 
within the soft parameters. When the flavour mixing is governed by MFV, it 
is possible for both the FV or 4B channels to be the dominant decay mode. 
Thus, it is important to study the signals of both possibilities at the LHC. 

We have investigated the limits placed by existing LHC searches at y/s — 
7 TeV on a light stop decaying entirely in the FV, 4B, or 3B channels. Our 
results suggest that these searches rule out a light stop decaying in these 
channels unless it is heavier than at least > 200 GeV. Even so, we 
emphasize that our results were obtained with an approximate simulation of 
the ATLAS and CMS detectors. We hope that a dedicated analysis will be 
performed by the collaborations to confirm and refine our findings. 

The exclusions we obtain also imply a significant bound on clcctroweak 
baryogenesis within the MSSM. To successfully create the observed baryon 
asymmetry, a very light, mostly unmixed stop is required with a mass less 
than that of the top quark. Our results suggest that such a state is ruled out 
by existing LHC analyses, at least if it decays promptly in the FV, 4B, or 
3B modes. Nevertheless, extensions of the MSSM that open up new decay 
channels might permit the stop to be as light as needed. 
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A MFV Conventions and Soft Parameters 

In this appendix we describe our detailed conventions for the Yukawa cou- 
phngs and the soft scalar parameters under the MFV hypothesis. We also 
describe how to match them to those used in the SLHA2 [HS] convention for 
specifying supersymmetric parameters. 

We define the Yukawa matrices appearing in the superpotential in the 
same way as Ref. (which follows Ref. [38]): 



(19) 



where we have written out the explicit SU (3)q x SU (2)„^ x SU (3)^^ indices. 
The relevant soft terms in the Lagrangian are 

+lJl{ai)%Q^-H^ - Dl{a\)\Q^.H, . (20) 



Both the superpotential of Eq. (19) and the soft terms of of Eq. (20) are 
invariant under the Gj = SU{3)q x SU{3)uh x SU{3)dii flavour group with 
the supermultiplets transforming as 



Q = (3,l,l), f/^ = (1,3,1), D= = (1,1,3), 
together with the Yukawa matrices transforming according to 

= (3,3,1), Yd = (3,1,3) . 



(21) 



(22) 



Suppose we make a field redefinition by a transformation (keeping the 
Yukawa matrices held fixed): 



Q -> VqQ, ^ U:U^ 



This modifies the effective Yukawa couplings appearing in Eq. (|19|) by 



M M ' V ' 



'^d ^ ^d 



(23) 



(24) 



In this picture, the MFV hypothesis for the soft masses is equivalent to 
the statement that the form of the soft terms given in Eqs. (|6}jT0|) remains 
unchanged under such a field redefinition, provided the old Yukawa couplings 
are replaced by the new ones {i.e. Yu Y^, Yd — )• in Eqs. (|6-10). 
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One of the two quark Yukawa matrices can be diagonalized by making 
Gf redefinitions. To diagonalize both simultaneously, we have to make a G/ 
transformation together with separate redefinitions for the two components 
of the Q = {uL^diY supermultiplet, 

ul -> Kml, cLl Vddi . (25) 

This gives 

^ f/KK = A„ (26) 

yI ^ UlYiVd = Xd, (27) 

where A„ = rriu/vu and A^ = nid/vd are diagonal with real, positive entries. 
In terms of the transformation matrices, the CKM matrix V is found to be 

V = VlVd . (28) 

This choice of field coordinates is called the super-CKM basis. Note that 
these combined transformations do change the form of the effective soft terms 
relative to the MFV forms of Eqs. ([6 10). 



To implement MFV soft terms, it is convenient to use G/ to choose a 
field basis where one of the two quark Yukawa matrices is diagonal before 
applying and Vd- The first and most popular choice is 

Yd = A,, Y^ = V^Xu , (29) 

where V is again the CKM quark mixing matrix. This is a nice choice because 
we only need to transform further hj Vu = V"^ (and Vd = 1) to diagonalize the 
quark masses. We use a second, similar choice in the discussion of squarks 
in Section [2} 

Yd = VXd, r„ = A, . (30) 

The additional transformation to diagonalize the quark masses is now Vd = 
VcKM (and K = !)• 

Non-diagonal scalar soft terms can be accommodated in the SLHA2 con- 
vention p9]. the down-type squark mass matrix in the super-CKM basis 
is 

^2^1 + + Dd^ Vdf^- fxmdtan(3 ^^^^ 

Vdfd - fJ'*md tan /3. m^^ + mj + 
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while the up-type squark matrix is 

^ f ^^l^^ + + Vufi - iirriu cot /3. ^^^^ 
y Vufu - iJ*mu cot p. mj^ + ml + D^^ 

Working in a basis where either Yu or is diagonal, these entries are 
related to the MFV soft masses using our conventions by 



yl^l^d = rh\ail + hV^XlV + b2Xl + . . .) (33) 
m}j = mlj. = m\a2l + hXl + ciX^VXlVXu) (34) 



= "^z)c = m^(a3l + &6-^d) (35) 

fu = aiVu = AXu^a^l + hVXlV^) (36) 

fa = a\Vd = AXd{a^l + hVXlV^). (37) 

These expressions apply for either Yu or Yd diagonal. 
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